Two-mode Bose-Einstein condensate in a high-frequency driving field that directly 

couples the two modes 
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A two-mode Bose-Einstein condensate coupled by a high-frequency modulation field is found 
to display rich features. An effective stationary Hamiltonian approach reveals the emergence of 
additional degenerate eigenstates as well as new topological structures of the spectrum. Possible 
applications, such as the suppression of nonlinear Landau-Zener tunneling, are discussed. An in- 
teresting phenomenon, which we call "deterministic symmetry-breaking trapping" associated with 
separatrix crossing, is also found in an adiabatic process. 
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I. INTRODUCTION 



Significant research efforts have been devoted to the 
nonhnear dynamics of interacting cold atoms, e.g., a 
Bose-Einstein condensate (BEG) in an external driving 
field. One main motivation is to understand how we can 
actively control the nonlinear dynamics and how the self- 
interaction of cold atoms can be used to simulate funda- 
mental models [l| . Examples include the control of the 
BEG self-trapping 0, i, Ij, effective turning-off of the 
self- interaction [5], controlled Mott-insulator transitions 
associated with a BEG in an optical lattice [6j], stabi- 
lization of bright BEG solitons by an oscillating mag- 
netic field tuned close to the Feshbach resonance [?!, as 
well as production of ultracoldjuolecules using stimu- 
lated Raman adiabatic passage [Oj. On a deeper level, 
BEG systems offer a useful tool to explore more aspects 
of many-body systems. In particular, the dynamics of 
a BEG in the large-particle-number limit is described 
by a mean-field nonlinear Schrodinger equation (Gross- 
Pitaeviskii equation). The resulting nonlinearity often 
challenges existing theories for linear systems. For ex- 
ample, the adiabatic following of a two-mode nonlinear 
system with an external field may necessarily break down 

Here we aim to examine how an adiabatic Landau- 
Zener (LZ) tunneling process of a BEG may be manip- 
ulated by an external driving field, thus extending an 
earlier study for linear systems [11]. As a second mo- 
tivation on a more fundamental level, we shall expose 
some nonlinear dynamics phenomena that do not exist in 
the mean-field dynamics of a BEG under field-free con- 
ditions. Specifically, we consider a two-mode BEG un- 
der a high-frequency field that directly couples the two 
modes (hence called "off-diagonal" driving below). A bi- 



ased static field is also considered for LZ processes. 

It is well-known that a high-frequency "diagonal mod- 
ulation" (e.g., high-frequency tilting of a double- well po- 
tential) can only re-scale the natural two-mode coupling 
strength [1, [ij [H, [l3|. However, our findings below 
for a different type of high-frequency modulation are dif- 
ferent. In particular, we show (i) that a high-frequency 
"off-diagonal" driving field can induce, in effect, addi- 
tional nonlinear terms in the mean-field equations of mo- 
tion, thus offering a control "knob" to tune different non- 
linear terms and simulate some systems not considered 
before; (ii) that different topological structures of the 
eigen-spectrum of the system can be generated and tuned 
by the driving field, also leading to additional degener- 
ate eigenstates. We then suggest to use a high-frequency 
driving field to realize the complete suppression of nonlin- 
ear LZ tunneling. In analyzing the adiabatic dynamics of 
a driven nonlinear system, we also find and qualitatively 
explain a phenomenon, called "deterministic symmetry- 
breaking trapping" associated with separatrix crossing. 



II. TWO-MODE SYSTEM UNDER 
OFF-DIAGONAL MODULATION 

The nonlinear two-mode system under off-diagonal 
modulation is described by 
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(1) 

where 7 denotes an external energy bias, |ap and |6p 
represent occupation probabilities for the two modes, c 
characterizes the nonlinear atom-atom interaction, and 
Aq denotes the static coupling between the two modes. 
We put h — 1 throughout this work. There are a num- 
ber of possibilities to experimentally realize this Hamilto- 
nian. For example, one may consider a BEG in a double- 
well potential, with the height of the potential barrier 
periodically modulated, or a BEG in an optical lattice 
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occupying two bands, with the well-depth of the optical 
lattice periodically modulated. In principle, these proce- 
dures should be achievable, considering previous experi- 
ments on two-mode EEC's [l^llBl- What might be even 
more feasible in realizing this two-mode system under off- 
diagonal modulation is to consider the internal states of 
a BEC, such as ^"^Rb [l3|, where there exist two internal 
states separated by a relatively large hyperfine energy. 
Then, the energy bias 7 can be effectively realized by the 
detuning of the coupling field from the resonance and the 
off-diagonal modulation may be achieved by modulating 
the intensity of the coupling field. Considering recent 
studies of two-mode nonlinear Schrodinger equations us- 
ing nonlinear optical waveguides (for example, see in Ref. 
Q), it might be also possible to realize our system in non- 
linear optics. 

Consider first the non-driven case, i.e. A = 0. Then 
H{t) reduces to the standard model of nonlinear LZ tun- 
neling (lo| . Therein the eigen-spectrum diagram as a 
function of 7 is known to display a loop structure at the 
tip of the lower (upper) level for c > Aq (c < — Aq). 
Such a loop structure, absent in linear systems, directly 
leads to a nonzero LZ transition probability even when 7 
changes adiabatically. As shown below, new system prop- 
erties emerge if the driving field is turned on. Without 
loss of generality we will restrict ourselves to the c > 
case, which requires an attractive interaction for bosons 
in a double-well potential or a repulsive interaction for 
bosons in two energy bands of an optical lattice. 

In the general case oi A with w 3> 7, c, Aq, it is 
found that |ap (= 1 — 1&|^) also oscillates at the frequency 
u. To expose possibly new physics hidden in the oscil- 
lations, another pair of wave function parameters (a', b') 
are found to be very useful, i.e., 

, _ « + -»^cos(^t) , Q " »^cosM) 
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b' = 



2 2 
Their equations of motion are given by 

'II ^ ^[7-cos(0) + c.cos2(0)(|6'p-|a'n 
-ic ■ sm{e) cos{e){a'*b' - a'6'*)]a' 
-hi[Ao + ■ sin(6') + c • sm^{e)ia'*b' 
-a'b'*) - ic ■ sm{d) cos(0)(|a'|2 - \b'\^)]b' 
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+ -[-7 • cos(0) - c • cos2(0)(|6'|2 - |a'|2) 
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(3) 



where 6 = ^cos(ti;t). Along with previous studies that 
focused on high-frequency driving fields [1, [13, [3 [13] j we 



consider now sufficiently large uj, such that the oscillation 
in is much faster than the natural time scale of the sys- 
tem as characterized by Aq, 7, and c (numerically, we find 
that the regime of w > 10 Aq, ui > 10c, and lj > IO70, 
where 70 is the initial value of I7I that is sufRciently large 
to ensure the LZ dynamics, can be safely regarded as a 
high-frequency regime; experimentally, a high-frequency 
driving field should not interfere with the two-mode de- 
scriptions). Then Eq. ([3]) can be significantly reduced by 
considering the averages of {a',b') over 2tt/uj. Speaking 
more rigorously, upon the large frequency condition, a 
zeroth-order approximation of a "1/a;" expansion can be 
used to yield 



.da' 
s— — 
dt 



db' 

i—r- 

dt 



= ^[7' + cz(l&? -!«?)]«' 

+ i[Ao + cy(a'*5'-a'6'*)]6' 
= i[Ao -cy(a'*6'-a'&'*)]a' 
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], and cy = c(sin (6*))^ 



] (Jo is the zeroth order Bessel function of 
the first kind). Evidently, these newly defined parameters 
reflect the action of the high-frequency driving field. We 
stress that the validity of this kind of high-frequency ap- 
proximation has been checked numerically and has been 
used in many situations. 

Equation (jH) no longer explicitly contains a time- 
dependent field. We can then define an effective static 
Hamiltonian Hcb that generates Eq. That is, 



7 + cz(|fop - |a|2) Ao + cy {a* b-ab*) 
Ao ~CY{a*b~ab*) -7 _ cz(|6|2 - |a|2) 



(5) 

where, for simplicity, we have replaced a' by a, b' by 
6, and so on. If we compare HcS with the original 
Hamiltonian in Eq. ([T]) for A = 0, one sees that the 
nonlinear parameter cz can be regarded as a rescaled 
parameter c, and the nonlinear term containing cy is 
new. In addition, the ratio of cz and cy is given by 
[l + Jo{2A/uj)]/[l — Jo{2A/uj)], easily adjustable by choos- 
ing different w and A. 

The effective Hamiltonian in Eq.® can be recognized 
as the one describing a single spin in a biaxial crystal 
field, with the cy (cz) term describing the anisotropy 
in the Y (Z) direction. Indeed, iJeff can also be writ- 



ten as Hsj 
Sz = 



= jSz + AqSx - czS 
, Sx = slbJ^ and Sy 



cySy, where 
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corresponding second-quantization Hamiltonian exactly 

describing the quantum system with TV bosons on the 

two modes is given by 
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FIG. 1: Upper panels: Level structures of the stationary ef- 
fective Hamiltonian H^s [see Eq. (O], as a function of 7. The 
dashed lines denote the new mean-field levels that are absent 
in a non-driven two-mode BEG. Symbols T, Dr, and Dl in- 
dicate how the involved levels are connected with the phase 
space structures shown in Fig. 2. Bottom panels: Parallel 
results in a fully quantum treatments for A'^ = 20. 



III. DETAILED RESULTS 

We now present in Fig. 1 the eigen-spectrum of Hcs 
as a function of 7. Evidently, the typical level structures 
(such as the loop structure) for a nonlinear LZ tunnel- 
ing model are also possessed by our system. On 
top of that, additional mean- field eigenstates (dashed 
lines) that are absent in a non-driven case also emerge 
through level bifurcations. The new eigenstates are di- 
rectly caused by the cy term induced by the driving field. 
In particular, if a loop structure exists and if cy < cz, 
then the additional level lies inside the loop, as shown 
in Fig. 1(b); and if cy > cz and cy > Aq, then level 
bifurcation takes place on the lowest branch and the ad- 
ditional level can be below the loop structure, as shown 
in Fig. 1(a). Figure 1(c) shows that the additional level 
may also exist in the absence of a loop structure. In the 
bottom panels of Fig. 1, we also show fully quantum me- 
chanical levels calculated from Eq. ^ for N = 20. The 
results confirm that the additional eigenstates we obtain 
on the mean-field level do have physical implications for 
fully quantum levels, even in the cases with a not very 
large N. 

Let us now examine Eq. ^ from a phase space per- 
spective, by mapping the mean-field trajectories of Eq. 
(ID) to that of a well-defined classical Hamiltonian system. 
The associated phase space can be defined in terms of s 
and (j), where (j) — cj)^ — (j)^ [11, [11], s = |6p — |ap, with 




FIG. 2: Phase space structures of He defined in Eq. (7), 
for Ao — 0.2; and cz = 0.4, cy = 0.6 (or equivalently, c = 1, 
A/uj = 1.42). "T" denotes an unstable fixed point with </> = tt. 
Dr and Dl denote two stable fixed points on the right and 
on the left. 

a = lale*"^" and b = |6|e*'^''. Using this pair of canonical 
variables, the involved classical Hamiltonian is: 

He - ^[-75- yS^-HAoVl-s^ cos(,^) 

_ ^^i^s')sm^m (7) 

The nonlinear eigenstates of Hcs now become fixed points 
in the phase space of He- Figure 2 displays phase space 
portraits of He for the parameters used in Fig. 1(a), for 
several values of 7 covering the regime of level bifurca- 
tion. In particular, the lower parts of Fig. 2(b), 2(c) and 
2(d) near (f> — n clearly show the splitting of one fixed 
point into three fixed points, thus associating the level 
bifurcation in Fig. 1 with the splitting of a fixed point. 
Because both the elliptic (stable) fixed points (marked 
by "Dij" and "Dl" in Fig. 2) yield the dashed line in 
Fig. 1(a), the additional level shown in Fig. 1 in fact de- 
notes two-fold degenerate eigenstates. By contrast, the 
hyperbolic (unstable) fixed point marked by "T" in Fig. 
2(d) yields the level right above the degenerate eigen- 
states [also marked by "T" in Fig. 1(a)]. In Fig. 2(f) 
and 2(g), the above-mentioned three fixed points start to 
merge back to one fixed point, in parallel with the level 
merging seen in Fig. 1(a) as 7 increases further. Exam- 
ining the phase space globally, it is also clear that the 
number of the fixed points and hence the number of the 
nonlinear eigenstates of H^s can vary from two to six, a 
clear sign that the nonlinear dynamics of a driven BEG 
can be very rich. 

It should also be noted that the above-mentioned two- 
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A right-moving trajectory 
crossing a swelling separatrix 



FIG. 3: A classical trajectory falling below an up-moving un- 
stable fixed point will evolve counter-clockwise around the as- 
sociated separatrix (solid line). The separatrix for a slightly 
larger 7 is represented by the dashed line. The shown trajec- 
tory is moving from left to right and will cross the swelling 
separatrix on the right and get trapped by the stable fixed 
point on the right. 

fold degeneracy occurs in a high-dimensional parameter 
space. In particular, for fixed nonlinear parameter c and 
fixed field parameters A and oj, the degeneracy can still 
occur in a two-parameter space of Aq and 7. This is in 
contrast to the well-studied non-driven model of a two- 
mode BEC where degeneracy occurs only along a line for 
fixed c, A and uj. 

So, how does the additional eigenstate shown in Fig. 
1(a) [cases in Fig. 1(b) and Fig. 1(c) are physically less 
appealing] affect the adiabatic dynamics? To answer this 
question we numerically solve Eq. ^ for the parameters 
used in Fig. 1(a), with 7 << and the initial state put on 
the lowest level. As 7 increases very slowly, the system's 
state is found to follow the non-degenerate lowest level 
up to the bifurcation point. When 7 increases beyond the 
"phase transition" point where the new two-fold degener- 
ate level emerges, the two-fold degenerate level becomes 
the lowest and the system is found to move along the 
new level. As 7 increases further, the two- fold degenerate 
level finally disappears and the system reaches the non- 
degenerate lowest level again, thus completing the LZ 
process. During the entire process, the system remains 
at the lowest level available and no transitions to any 
upper levels are found. Hence, the new two-fold degen- 
erate level induced by the driving field offers a means to 
circumvent the loop structure and hence totally suppress 
the nonlinear LZ transition that is doomed to happen if 
the two-fold degenerate state were not there. To con- 
nect this observation of totally suppressed LZ tunneling 
in the (a', 6') representation [see Eq. (2)] with the direct 
observable (a, b) in experiments, note that (i) initially if 
A = then a = a' and b = b', and (ii) if A is switched 
on slowly enough as compared with oj but fast enough as 
compared with the change rate of {a',b') (characterized 
by 7, Aq and c), then the initial values of (a, b) are passed 
to (a',6')- 

One more aspect of the above nonlinear LZ process 



remains to be examined. Because the dashed line in 
Fig. 1(a) denotes a two- fold degenerate eigenstate, we 
should study which state the system will reside in when it 
slowly pass the level bifurcation point with an increasing 
7. Since the phase space structure shown in Fig. 2 always 
possesses a mirror symmetry with respect to 4> = n, one 
may intuitively expect that during the LZ process the 
system is trapped by either of the two stable points Dji 
or Dl in a random fashion, with equal probability. How- 
ever, we find that this picture is incorrect here. Instead, 
the system is found to be deterministically trapped by 
Dfi [see Fig. 2(d)]. Physically, this deterministic trap- 
ping means that the relative phase between a' and b' is 
not random during the LZ process, i.e., it is robust to 
small fluctuations in the initial state. A careful analysis 
enables us to explain this intriguing observation qualita- 
tively. As 7 increases, the fixed point with (j) = tt [see the 
lower part of Fig. 2(c) and Fig. 2(d)] moves upwards in 
the phase space. When this fixed point becomes unstable 
[denoted "T" in Fig. 2(d)], the adiabatic following must 
break down and hence the actual trajectory will find it- 
self slightly below the up-moving fixed point. As such, 
the trajectory starts to slowly move counter-clockwise 
around a separatrix, or from left to right, as illustrated 
in Fig. 3. At the same time, because 7 is increasing, 
the separatrix deforms and swells, and as a result the 
trajectory necessarily crosses the separatrix on the right 
and hence gets trapped by Dn. Indeed, if we reverse the 
adiabatic process, i.e., passing the bifurcation point with 
a decreasing 7, then one can predict that the system will 
be trapped by D^, a prediction confirmed numerically. 
This counter-intuitive "deterministic symmetry-breaking 
trapping" is complementary to a well-known separatrix- 
crossing-induced phenomenon, i.e., "quasi-random" trap- 
ping in classical mechanics (which has been systemati- 
cally applied to a few BEC systems @). Encouraged by 
the finding here and in efforts to confirm the generality 
of our finding, we also studied the adiabatic following 
dynamics of a modified rotating pendulum system whose 
fixed point moves with an external parameter. Analogous 
results are also found in this pendulum case. Hence, it 
should be of some interest to carry out an experimental 
BEC study of the observed symmetry-breaking separa- 
trix crossing here. Such kind of experiments might also 
offer new insights into the validity of the mean-field de- 
scription of a BEC. 



IV. CONCLUSION 

To conclude, we have theoretically examined the dy- 
namics of a two-mode BEC driven by a high-frequency 
driving field that directly couples the two modes. Based 
on our results here we expect rich phenomena in general 
for a multi-mode BEC under a high-frequency driving 
field. Though our results are purely theoretical, it is our 
hope that the results here will stimulate future experi- 
ments on the dynamics of a BEC in high-frequency driv- 
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ing fields and on the control of nonlinear LZ tunneling 
dynamics. 
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